Realization of a quantum walk with one and two trapped ions 
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We experimentally demonstrate a quantum walk on a line in phase space using one and two 
trapped ions. A walk with up to 23 steps is realized by subjecting an ion to state-dependent 
displacement operations interleaved with quantum coin tossing operations. To analyze the ion's 
motional state after each step we apply a technique that directly maps the probability density 
distribution onto the ion's internal state. The measured probability distributions and the position's 
second moment clearly show the non-classical character of the quantum walk. To further highlight 
the difference between the classical (random) and the quantum walk, we demonstrate the reversibility 
of the latter. Finally, we extend the quantum walk by using two ions, giving the walker the additional 
possibility to stay instead of taking a step. 



The Galton board [l[ is a mechanical device in which a 
falling ball encounters a triangular lattice of pins stuck in 
a board that repeatedly scatter the ball to the left or right 
in a random way. Originally conceived for illustrating 
the emergence of normal probability distributions, it can 
also be considered as an apparatus for carrying out a 
random walk on a line [1, Q , a notion that had not been 
introduced into the scientific literature at that time. 

Since then, random walks have become an ubiquitous 
concept in physics and computer science. The quantum 
walk [3, H| is the quantum analogue of a random walk. 
In its discrete one-dimensional version, a spin- 7; quan- 
tum particle initially described by a wave packet cen- 
tered at position xq undergoes a one-dimensional motion 
governed by the particle's internal state. The particle is 
state-dependently displaced by a step of length d by the 
action of the unitary operator Ud = exp(—icrjpd) where 
a j is a spin projection operator and p the momentum op- 
erator (see Fig. 1). This operation is followed by another 




unitary operation Ui = exp(— i\(Jk) with Tr(ajak) = 
scrambling the particle's internal state. After N itera- 
tions of this elementary step, the particle's initial wave 
function l^o) has evolved into 
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After TV steps, for a wave packet initially localized at xq = 
0, the wave packet is spread out over a distance 2Nd. 
Moreover, due to quantum interference of different paths, 
the spatial probability distribution strongly differs from 
the classical case. While for the classical random walk a 
binomial probability distribution results with (a; 2 ) cx N, 
the distribution for the quantum walk is peaked towards 
the outer edge and has a second moment growing like 
(x 2 ) oc iV 2 . 

There have been a number of proposals discussing 
experimental realizations of one-dimensional quantum 
walks in systems like atoms in optical lattices @ , trapped 



FIG. 1: Quantum walk in phase space. In each step of 
the walk, a state-dependent displacement operation splits the 
wave function in phase space into two parts followed by a coin 
tossing operation that coherently scrambles the internal state 
of the ion. These operations are repeated N times. To mea- 
sure marginal distributions in phase space, a probe pulse is 
applied that state-dependently displaces the wave function in 
phase space in a direction orthogonal to the one to be mea- 
sured. 



ions 0, or cavity QED Q. Recently, experimental real- 
izations with atoms in an optical lattice a trapped 
ion [l(| , and photons 11 1 have been reported. 

For the case of trapped ions, different techniques for 
analyzing the quantum walk have been discussed @, El] 
and a proof-of-principle experimental realization was re- 
ported recently [l(| for a limited number of steps. In this 
paper, we demonstrate a discrete quantum walk with up 
to 23 steps using a single trapped ion and analyze it by 
a measurement technique that directly reconstructs the 
ion's probability density along a line in phase space. 

In the experiment, a single 40 Ca + ion is suspended in 
a linear Paul trap [13f with radial and axial trap fre- 
quencies of uj r « (2ir) 3 MHz and w ax = (2tt) 1.356 
MHz, respectively. Doppler cooling, resolved sideband 
cooling of the axial mode and optical pumping prepare 
the ion in the ground state of motion and the internal 

21] . A narrow linewidth 



state |Si/ 2 , m = 1/2) = |- 
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laser at 729 nm coherently couples the states \—) z and 
I-D5/21 m = 3/2) = \+) z . State detection is done via flu- 
orescence detection on the S-y/2 <-> Px/2 transition (l3| . 

A general state-dependent displacement Hamiltonian 
is implemented using a Dichromatic light field at 729 nm 
that is resonant with both the blue and red axial sideband 
of the \—) z <-> \+)z transition. In the Lamb-Dicke regime, 
the resulting Hamiltonian, which is the sum of a Jaynes- 
Cummings and an anti-Jaynes-Cummings Hamiltonian, 
is given by 



Hd =Hr)Cl {{u x cos 4>+ — &y sin 4> + ) 

<8> ((a + a' ) cos 0_ + i(a' — a) sin </>-)) • (2) 

Here, 77 = 0.06 is the Lamb-Dicke parameter, ft the Rabi 
frequency and 2</> + = 4> r + (f>^ and 2</>_ = 0b — 4> T are the 
sum and the difference, of the phases of the light fields 
tuned to the red and blue sideband. 

To perform a symmetric quantum walk the ion is pre- 
pared in the state \+) y = (\+) z + i\~)z)/V2 by a n/2- 
pulse on the carrier transition. Applying the bichromatic 
light field with </>_ = ir/2 and <f> + = realizes the Hamil- 
tonian Hd = 2i]flA x a x p with the momentum operator 
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Hamiltonian for a duration r generates the propagator 
Ud with step size d = 2tjD,tA x . 

Under the action of Hd, the ion's wave packet coher- 
ently splits in phase space along the x-axis. The two 
emerging wave packets "01 1 V4™^ are associated with 
the internal states |±)a;. The length and the intensity of 
the pulse determine the width of the splitting. In our ex- 
periments we use a pulse of 40 ixs with a Rabi frequency of 
51 = (2tt) 68 kHz to achieve a step size of d — 2A X . This 
step size makes the two resulting motional wave packets 
nearly orthogonal, KVq 1^4 )| 2 ~ 0.02, but still allows 
for a large number of steps in phase space. Next, we 
perform a 7r/2-pulse acting on the carrier transition as a 
symmetric coin flip. This pulse creates an equal super- 
position of a x eigenstates for both wave packets. These 
two pulses are repeated according to the number of steps 
to be carried out. 

To measure the probability distribution along a line 
in phase space, we create two displaced copies of the 
state that are subsequently interfered. For this, we 
use of another state-dependent displacement operation 

measurement of a z fol- 
lowing the application of U p is equivalent to measuring 
the observable 



O(k) = Ul<j z U p 



cos(kx)a z + sin(kx)a y 



(3) 



with the usual position operator x = (aS + a)A x on the 
initial state. The propagator U p is obtained by setting (f> + 
and 4>- in Hd to 0. Here, k = 2rfX v tj A x is proportional 
to the interaction time t. If the ion's internal state is 
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FIG. 2: (a) Measurement of Fourier components (cos(kx)) 
and (sin(fca;)} for a seven-step quantum walk. The data are 
obtained by varying the duration of the probe pulse for the 
ion prepared in the internal state |+) z (left) or |+) y (right) 
after completing the walk. The probability distribution is ob- 
tained by Fourier transforming a fit to the data (solid line), 
(b) Reconstruction of the symmetric part of the probability 
distribution (5(x — x)) for up to 13 steps in the quantum walk. 
The blue dashed curve is a numerical calculation for the ex- 
pected distribution within the Lamb-Dicke regime. The blue 
solid curve takes into account corrections to the Lamb-Dicke 
regime. In step 7, the dotted curve represents the full recon- 
struction using also the (sm(fcr)) shown in (a).(c) Probability 
distribution of a five-step quantum walk after application of 
five additional steps which invert the walk and bring it back 
to the ground state. 



|+)z, we have (0(k)) = (cos(fci)) and for |+)j,, we have 
(0(k)) = (s'm(kx)). A Fourier transformation of these 
measurements yields the probability density (S(x — x)} 
in position space which for a pure state |^) amounts to 

|\I/(:e)| 2 . Furthermore, we have that ^s(0(k)) 

(x 2 <7 z ) (l6j |. For eigenstates of a 
the expectation value (0(k)) thus gives the width of the 
probability distribution w x . 

The quantum walk entangles internal and motional 
degrees of freedom. Its analysis, however, requires the 
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FIG. 3: (Color online) (a) Width w x of the probability dis- 
tribution in units of ground state size A^ as a function of the 
number of steps for a quantum (■) walk. The solid curve 
represents a full numerical simulation of the quantum walk as 
realized in the experiment. The width of the x-distribution 
for a classical random walk (•) increases more slowly and is 
described (solid red line) by eq. (J4j . The data points (♦) show 
the measured width w p of the marginal distribution along the 
p-direction with A p = h/2A x . (b) Average number of vibra- 
tional quanta after N steps in the quantum walk measured 
by driving oscillations on the carrier transition. The solid 
line is based on a full simulation, the dashed line assumes the 
validity of the Lamb-Dicke approximation. 



preparation of pure internal states like \+) z or \+) y - 
Therefore, we recombine all internal state populations 
in \—}z before the measurement. To this end, the popu- 
lation in is transferred to |— } z after transferring the 
population in \—) x to the auxiliary state I-D5/2, m = 5/2). 
A laser pulse at 854 nm excites the population from 
l^5/2i m = 5/2) to I-P3/2,™ = 3/2) from where it spon- 
taneously decays to |— ) z . The efficiency of this pumping 
process is > 99 %, limited by a small branching ratio 
to the _D 3 / 2 -state. Only after the recombination step, 
we prepare the internal state required for measuring the 
even or odd Fourier components of (|5|) . Due to the small 
Lamb-Dickc parameter, the probability of changing the 
motional state of the ion during the pumping steps is 
small and hardly affects measurements of observables in 
position space at all. 

In the experiment we set fl p = (2tt) 26 kHz and mea- 
sure (a z ) for probe times between and 300 fjs in order 
to reconstruct the probability distribution (S(x — x)} for 
different numbers of steps N. Since the walk is sym- 
metric, it is in principle sufficient to measure only the 
even components of ([5]). For a seven-step walk, the mea- 
sured odd and even Fourier components are displayed in 
Fig.[2[a). Panel (b) shows the reconstructed probability 
distribution (S(x — x)) based on the even components for 
up to 13 steps. The uneven terms were checked to be 
close to zero for each number of steps N. The dashed 
lines in the plots are numerical simulations based on the 
Lamb-Dicke approximation. These lines deviate from the 
reconstructed distribution for N > 7 due to higher order 
terms in 77 that are not taken into account in eq. (|9|). 
The solid lines are based on a numerical simulation us- 
ing all orders. A similar difficulty occurs in the mea- 
surement of observables based on eq. ([S]). For this rea- 
son, the reconstruction is not accomplished by a direct 



Fourier transformation of the data. Instead, we apply 
a constrained least-square fit based on convex optimiza- 
tion 17[ capable of handling higher-order corrections (see 
the EPAPS document for more information on the re- 
construction process). To get smoother distributions ad- 
ditional constraints were invoked by the reconstruction 
algorithm. A physical constraint is given by the max- 
imal kinetic energy a one-dimensional wave packet can 
have. An estimate for the kinetic energy can be deter- 
mined by measuring the momentum distribution in the 
same way as the position distribution. By changing </>_ 
to 7r/2 in the probe pulse, the operator x appearing in 
([5|) is replaced by an operator oc p. These measurements 
(see Fig. [3](a)) indicate that the momentum distribution 
is not seriously affected during the walk, as expected for 
a pure displacement along the x-axis. 

A striking difference between classical and quantum 
walks is the reversibility of the latter. In the experiment 
we reversed a quantum walk after five steps. This was 
done by switching the phase of the following five dis- 
placement and coin flipping pulses by 7r. In this way the 
quantum walk is exactly reversed and the ion returns to 
the ground state. The corresponding reconstructed prob- 
ability distribution shown in Fig. [2jc) closely resembles 
the one of the initial state and demonstrates once more 
the coherence of the quantum walk. 

To further highlight the differences between quantum 
and classical walks we also realized a classical walk by 
randomizing the phase between each step (while keeping 
the coin flip-displacement operator pair coherent for each 
individual step). The phase for each step was generated 
by a random noise generator. This mimics a completely 
mixed ensemble of measurement outcomes that behaves 
classically. A good way of quantifying the difference be- 
tween the quantum and classical walks is by measuring 
the average width of the probability distributions. For a 
classical walk with a step size d — sA x we have 



w x = A a 



2s N 



1. 



(4) 



where the second term takes into account the initial 
width A;,, of the probability distribution. By contrast, 
for a quantum walk the width goes as w x ~ N for high 
N. To measure w x for the random walk, the curvature of 
(<j z ) at short probe-time was analyzed. Quadratic fitting 
gives direct access to the width w x . For the quantum 
walk, w x was obtained from the measured probability 
distributions. In Fig. |3|a) the results of these procedures 
can be seen seen for both a quantum and a classical walk. 

To avoid problems in the measurement of the motional 
state due to leaving the Lamb-Dicke limit for large num- 
bers of steps, we implemented a method suggested in 12 1 . 
Outside the Lamb-Dicke regime the coupling strength 
Q n ,n on the carrier depends on the phonon number n 
as fl n , n = QoLniv 2 )- Here, L n (rj 2 ) is the n-th order 
Laguerre polynomial. The mean phonon number (n) is 
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FIG. 4: Reconstructed probability distribution {S(x — x)) for 
a two-ion quantum walk with up to 5 steps with a step size 
of 4A*. 



computation [19( and in finding search quantum algo- 
rithms that outperform their classical counterparts [20(. 
As such the experimental implementation of the quan- 
tum walk serves as an important benchmark and points 
the way to further experiments. For instance, the imple- 
mentation of a quantum walk with two ions opens up the 
interesting possibility to introduce entanglement [18[ and 
more advanced walks. 



Appendix: Reconstruction of the probability density 

For the reconstruction of the probability density 
p{x) = (S(x — x)) of the motional quantum state p m , 
we determine the expectation value of the observable 

O(k) = Up<j z U p = cos(kx)a z + sin(kx)a y (5) 



determined by a constrained least-square fit of the car- 
rier Rabi flops with the number state distribution as a fit 
parameter. In Fig.^b), the resulting average vibrational 
quantum numbers are shown. As expected for the quan- 
tum walk, we observe a quadratic dependence (n) cx iV 2 
on the number of steps. 

Finally, we extend the quantum walk concept by 
adding a second ion to the system [18| . In the two ion 
quantum walk we make use of the center-of-mass mode. 
To account for the second ion, all Pauli matrices <x, in 



cq. |T]) are replaced by a\ + <r\ . This changes the coin 
from two sided to four sided, with three possible opera- 
tions. The "side" belonging to the state | + +) x (| ) x ) 

corresponds to a step to the right (left) while the sides 
belonging to the states | H — ) x and | — \-) x correspond 
to no step at all. The ions are prepared in the state 
I + with a 7r/2-pulsc leading to a symmetric walk. 
For the two ion quantum walk all pulses arc applied to 
both ions simultaneously. The probability distribution of 
the center of mass mode is obtained in the same way as 
for a single ion. The results for a walk of up to 5 steps are 
shown in Fig. 2) Again, the distribution deviates strongly 
from the classical version and shows a faster spreading. 

In summary, we have implemented a quantum walk us- 
ing trapped ions. An experimental technique was imple- 
mented to determine the probability distribution along 
a line in phase space. This method might have further 
applications in quantum optics experiments or quantum 
simulations [15j. We have highlighted the difference be- 
tween a classical and a quantum walk and demonstrated 
the reversibility of the latter. The current limitation in 
number of steps is given by instabilities in the trap fre- 
quency leading to decoherence and by the change in the 
coupling strength due to high phonon numbers. Quan- 
tum walks are of importance as a primitive for quantum 



by applying the unitary U p = exp(— ikxa x /2) to the 
state p = |^ r )(^ r | (8 Pm and measuring the operator 
a z . The right-hand-side of eq. (|5|) is obtained by us- 
ing the equality exp(i9a x ) = cos 9 1 + is'm9a x and 
<7,(7j = CijkO'k for i 7^ j. In this way, we determine 
(cos(fci)) by choosing \^>) = \+) z , and (sin(fcr)) by 
|^) = \+) y 14] . In principle, a Fourier transformation of 
f(k) = (cos(fci)) + i(sin(fcr)) is sufficient for obtaining 
the density (S(x — x)) in position space. However, with a 
finite number of experiments, the expectation values can 
only be determined for a discrete number of /c-values and 
these measurements do not yield the exact expectation 
values but rather estimates of them. As a consequence, 
a reconstruction of the density based on Fourier trans- 
formation gives unphysical probability densities that are 
not non-negative everywhere. To overcome this prob- 
lem, we reconstruct p(x) by a constrained least-square 
optimization based on convex optimization 17j . We dis- 
cretize the position space by using a suitable set of points 
Xi and search among the probability distributions with 
p{xi) > for all Xi the distribution that minimizes 



S = ^ (^2p{x l )cos(kx l ) - C^j 



(6) 



where Ck and Sk are the experimentally determined es- 
timates of (cos(fcr)) and(sin(fc£)), respectively. 

In our reconstruction, we use another physical con- 
straint based on a measurement of the kinetic energy 
1 ^—^ in the following way. For a wave function -0 = 



\2m> 

A(x)e 



4(x) 



a lower bound on the kinetic energy is given 
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by 



= dx((A(x)^(x)f+A'(xf) 
2m 2m 



f.2 roo 
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~ 2mj_ 00 v ' 



(7) 



where differentiation with respect to x is indicated by 
primes. For p(x) = \ip(x)\ 2 , we then have because of 
A(x) = p{x)? and A'(x) = \p'{x)p{x)~ ' that 



2m 8m J _ OQ p(x) 



(8) 



This constraint is also valid for mixed quantum states. 
In our optimization algorithm, cq. ((5J) is a convex con- 
straint that excludes distributions p(x) having excessive 
energies. It requires a measurement of (p 2 ) which is ob- 
tained by setting </>_ = n/2 in the bichromatic Hamilto- 
nian generating the unitary U p in eq.® and calculating 
d 2 1 'dk 2 (0(k)) . Adding the constraint works particularly 
well for the states produced by the quantum walk. These 
states ideally do not have any phase gradients 4>'(x) in 
which case inequality (JSJ turns into an equality. 

Moreover, the optimization algorithm can also handle 
to some extent problems related to the validity of the 
Lamb-Dicke approximation. In this approximation, the 
bichromatic laser-ion interaction, which is used for the 
operation U p in the reconstruction measurement, is de- 
scribed by the Hamiltonian 



Hd = hrjVta x ® (a + a?). 



(9) 



This Hamiltonian is strictly valid only for rj — > because 
it is based on a Taylor expansion g i v( a + ai ) = J + i-qia + 
a') + 0(t] 2 ) of the atom-light interactions that neglects 
terms in r\ of order two or higher. If resonant terms up 
to third order are taken into account, the Hamiltonian 
becomes 



77 

(a + a f ) - ((a + a f )n + h(a + a f ) + l) 



Since this Hamiltonian no longer commutes with x, it 
cannot be used instead of eq. (fTOj) for the reconstruction 
procedure described above. On the other hand, an anal- 
ysis based on eq. ([§]) yields wrong results for quantum 
walks for large number of steps where the created states 
no longer fulfil the Lamb-Dicke criterion. Fortunately, for 
these states, their potential energy ex ((a + a.t) 2 ) is much 



larger than their kinetic energy oc ((i(a — a^)) 2 ) which 
provides the justification for replacing n by j(a + a*) 2 in 
cq. (fll)]) . Using this assumption, the Hamiltonian 



Hjy = fvqVL a x 



( a + a t)(l_ZL(( a + a t)2 + 1 V ) 



(11) 



becomes again diagonal in the real space basis. The re- 
constructed probability densities (5(x — x)} shown in the 
paper are based on this Hamiltonian. 
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